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Abstract 

In this study, the temperature variations of the equilibrium and the non- 
equilibrium antiferromagnetic and ferromagnetic susceptibilities of a metamagnetic 
system are examined near the critical point. The kinetic equations describing the 
time dependencies of the total and staggered magnetizations are derived by utilizing 
linear response theory. In order to obtain dynamic magnetic relaxation behavior of 
the system, the stationary solutions of the kinetic equations in existence of sinusoidal 
staggered and physical external magnetic fields are performed. In addition, the static 
and dynamical mean field critical exponents are calculated in order to formulate 
the critical behavior of antiferromagnetic and ferromagnetic magnetic response of 
a metamagnetic system. Finally, a comparison of the findings of this study with 
previous theoretical and experimental studies is represented and it is shown that a 
good agreement is found with our results. 
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1 Introduction 



Investigation of the AC or dynamic susceptibilities are the most preferred 
techniques to study Magnetic relaxation (MR) which is related to magnetic 
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hysteresis and present in all stages of the development of ordered phase is an 
important approach to probe magnetic systems. They are obtained from the 
dynamic response of the system to time-dependent magnetic field and now 
commonly used to investigate the magnetic properties of high-T c systems [1] , 
cobalt-based alloys [2] , nanoparticles [3] , spin glasses [I] , molecule-based mag- 
nets [5], and magnetic fluids [6j. Theoretical investigation of dynamic magnetic 
response of the Ising models has been an subject of interest for quite a long 
time: In 1966, Barry has studied spin-1/2 Ising model by a method combin- 
ing statistical theory of phase transitions and irreversible thermodynamics [7] . 
Using the same method, Barry and Harrington has focused on the theory of 
relaxation phenomena in an antiferromagnet [8]. In addition, they obtained 
the temperature and frequency dependencies of the magnetic absorption and 
dispersion factor in the neighborhood of the critical temperature. On the other 
hand, Suzuki and Kubo has obtained the time dependent susceptibility of the 
kinetic Ising model [9]. Acharyya and Chakrabarti presented the real and 
imaginary parts of magnetic susceptibility near the orderdisorder transition 
point of a spin-1/2 Ising system in the presence of a periodically varying ex- 
ternal field by using Monte-Carlo simulations [TU]. On the other hand, The 
dynamic magnetic response of the materials and the development of meth- 
ods for its modification are important for their potential applications: Cores 
made of cobalt-based alloys in low signal detectors of gravitational physics 
contribute as a noise source with a spectral density proportional to the ac 
susceptibility of the alloy [2]; MR effects in nanocomposite particles are used 
in design of magneto-optical devices the superposition principle for the 
imaginary part of complex magnetic susceptibility of composite magnetic flu- 
ids is crucial for the design of absorbers and microwave attenuators which are 
based on ferromagnetic resonance absorption of the electromagnetic field [6]. 

Recently, Erdem investigated the magnetic relaxation in a spin-1 Ising 
model near the second-order phase transition point. In this study, time deriva- 
tives of the dipolar and quadrupolar order parameters are treated as fluxes 
conjugate to their appropriate generalized forces in the sense of irreversible 
thermodynamics In addition, Erdem has obtained the frequency depen- 
dence of the complex susceptibility for the same system [TJ]. Though the 
MR has been subject of many theoretical and experimental investigations 
mentioned above, there has been no study to investigate the behavior of ac 
susceptibility in the metamagnetic Ising system. 

Metamagnets, systems in which antiferromagnetic and ferromagnetic in- 
teractions exist simultaneously, are of great interest because it is possible to 
induce novel kinds of critical behavior by forcing competition between these 
couplings, in particular by applying a magnetic field [13] . The metamagnetic 
model that has in-plane ferromagnetic coupling and antiferromagnetic cou- 
pling between adjacent planes, and the next nearest neighbor (nnn) spin-1/2 
Ising model with antiferromagnetic nearest-neighbour (nn) and ferromagnetic 
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nnn interactions are the theoretical Hamiltonian models showing a similar kind 
of behavior. These models have been investigated by Monte-Carlo simulations 
[H], as well as High-temperature series expansion calculations [15] . 



FeBr 2 and FeCli are typical Ising type metamagnets [13]. In these struc- 
tures, in the antiferromagnetic phase when the iron ions in the triangular 
layers order ferromagnetically, a layer with a negative sign follows a layer with 
a positive sign. Due to this fact, the external field acts differently on the oppo- 
sitely oriented layers which leads to different ordered states and associates a 
sequence of phase transitions as a function of these two interaction strengths. 
A Monte-Carlo simulation has been performed on a realistic model of FeC/2 
under an external magnetic field [16J, in addition this typical metamagnet 
has also been treated by a high-density expansion method on a two-sublattice 
collinear Heisenberg Ising model with three- and four-ion anisotropy [T7fT8] . 
Dynamic properties of these systems have been investigated by using three- 
spin flip dynamics [T9] , Kawasaki dynamics [20], Glauber dynamics [21], and 
dynamic Monte Carlo renormalization group method [22] . In our recent works 
we have calculated the kinetic phase diagrams of the system under an oscillat- 
ing field [23.24] . Moreover, we have presented an investigation of the relaxation 
dynamics of iron group dihalides and studied the field and temperature de- 
pendence of the relaxation times via the phenomenological kinetic coefficients 
[25] . Recently Gulpinar and Karaaslan performed an investigation of the re- 
laxation times with taking in account the interference between the relaxation 
processes of antiferromagnetic and ferromagnetic order parameters [26] . 



The purpose of this paper is to study the dynamical magnetic response 
properties of the spin-1/2 Ising model in the presence of oscillating external 
magnetic field. On the other hand, to the best our knowledge there is no 
study in the literature which represents the mean field calculations for the 
static staggered and direct susceptibilities of a metamagnetic Ising system. 
Since then we have given the derivation of the static magnetic response func- 
tions and their temperature variations near the critical point in Sec. 2. Then, 
in Sec. 3 we derived the temperature and frequency dependent dynamic stag- 
gered and direct susceptibilities by a similar method previously used to study 
relaxation dynamics and sound propagation [27] and investigate the behaviors 
of magnetic dispersion and absorption factors near the second order phase 
transition temperature. The method used in this paper provides information 
about the dynamical critical properties based on the phenomenological ki- 
netic coefficients of the phenomenological rate equations, which governs the 
magnetization relaxation, is due to an "ad hoc" spin-lattice coupling. 
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2 Derivation of Static Staggered Magnetic and Magnetic Suscep- 
tibilities of Spin-1/2 Metamagnetic Ising Model 



In order to obtain static staggered magnetic susceptibility one should 
introduce a staggered external field B s to the system [28] whereas magnetic 
susceptibility is the response of the total magnetization to a physical exter- 
nal field B. Since then one should add two different magnetic fields to the 
Hamiltonian of the spin-1/2 Metamagnetic Ising Model: 



H = - J E s i S i + E Si Sj - fi B J2 S { - fiB s ( E St - E s, ) 

nn nnn i=l \ 1 2 / 



(1) 



J, J ,n,B and B s are ferromagnetic and antiferromagnetic exchange interac- 
tion constants, spin magnetic moment, external magnetic field and external 
staggered field respectively. By making use of mean field approximation free 
energy of the system can be obtained as following: 

G(m t , m s , T, B, B s ) = \ JN{m 2 t - m 2 s ) - \j' N{m 2 s + m\) 

-\k B TN(A ln(2) - (1 + m s + m t ) ln(l + m s + mt) 
— (1 — m s — m t ) ln(l — m s — m t ) — (1 + m t — m s ) 
ln(l + m t — m s ) — (1 — m t + m s ) ln(l — m t + m s )) 
-Ngfi B Bm t - Ngfi B m s B s , 

(2) 

where k B , g, ^b-, N are the Boltzmann's constant, the spin factor, the Bohr 
magneton and the total number of metamagnetic Ising spins respectively. For 
a simple cubic lattice in which intralayer interactions are ferromagnetic and 
interlayer interactions are antiferromagnetic z = 2 and z — 4. The equilib- 
rium conditions, ( = 0, = ) result in the following transcendental 
\dm t om s J 

equations: 



m t 



m. 



sinh(2^; 



cosh(2^) + cosh(2^H)' 

(3) 

sinh(2^±f^) 



cosh(2^) + cosh(2^H; 

Here, / = J — J' and t = J + J'. These equations may be solved without 
difficulty by an iterative procedure, i.e. Newton Raphson method. The equi- 
librium solutions should correspond to extremum of the free energy so that 
one has to determine the solution that minimizes G [28\. Since the solution 
of these equations are discussed in Ref. [2S] extensively, we shall only give 
a brief summary here as follows: Topology of the metamagnetic Ising model 
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phase diagram depends on the value of the ratio of the exchange interactions 
(77 = jj): (i)If 77 > 0.6, the transitions between the anti-ferromagnetic and 
paramagnetic phases are of first order at low temperatures and strong fields 
while it is of second order at higher temperatures. The two types of transi- 
tions are connected by a tricritical point, (ii) For < rj < 0.6, the tricritical 
point decomposes into a critical end point (CEP) and a double cricital end 
point (DCP) with a line of first order transitions in between, separating two 
anti-ferromagnetic phases [28125], see Fig.l(b)-(c) in Ref.[2S]. By definition 
staggered magnetic susceptibility is, 

Xstt - = lim (B = B ) (4) 

and total (direct) magnetic susceptibility can be expressed as, 

Xttr = lim ^r 1 (B = B ). (5) 



If one uses the equations of state given in Eq.(|3j), the static staggered magnetic 
and magnetic susceptibilities can be found as, 



lim 



a 2\C\ — aii c 2 



B s ^o \a 22 a n - ^21^12 



Yt * *■ = lim 

/V L static tj . 



(6) 

622^1 - 612^2 \ 



B s ^0 \ 611622 - 612621 



Where, a n ,ai2, a 2 i, a 2 2, ci, c 2 , and 6n, 612, 6 2 i, 622, c?i, d 2 are given in Appendix- 
A. 

Fig.l(a)-(b) represents the critical behavior of static staggered magnetic 
susceptibility (xs) an d direct magnetic susceptibility. The arrows illustrate 
the phase transition temperature Tq{B) for B = 1.0. The staggered magnetic 
susceptibility (xs) increases rapidly with increasing temperature and diverges 
as the temperature approaches to the second order phase transition point on 
either side, as seen in Fig. 1(a). On the other hand, the magnetic suscepti- 
bility (xt) also increases rapidly when the temperature is raised but makes 
finite jump discontinuity at the second order phase transition temperature, 
which is illustrated in Fig. 1(b). These findings are in accordance with the re- 
sults given in Refs. pOlBT] obtained within effective field approximation for an 
diluted metamagnetic Ising Model. We should also note that Barry and Har- 
rington have given an detailed analysis for the static total susceptibility of an 
antiferromagnetic Ising system in the Bethe-Takagi approximation which has 
the same result of the constant-coupling approximation of Kasteleyn and Van 
Kranendock applied to Ising antiferromagnetism [32]. In both studies an finite 
jump discontinuity is observed at the critical points which is in accordance 
with our results. On the other hand, a infinite-slope singularity is found by 
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exact series expansions [33]. One might also emphasize that the mean field 
approximation only barely fails to exhibit a maximum in the static direct sus- 
ceptibility above the T/y, that implies a need for a method with takes into spin 
correlations [31] . 



3 Derivation of Non-equilibrium Staggered Magnetic and Mag- 
netic Susceptibilities of Spin-1/2 Metamagnetic Ising Model 

In order to study the relaxation processes in the anhydrous iron 
group dihalides, within the spin-| metamagnetic Ising model, one assumes a 
small deviation in the value of the applied external magnetic field. This fact 
removes the system slightly from its equilibrium state, and one can investigate 
how rapidly the metamagnetic system relaxes back to its equilibrium state. 
On the other hand, it is well known that the metamagnetic - paramagnetic 
phase transition lines occur at places which are away from the B = axis 
[28] . Consequently, the magnetic Gibbs free energy production (AG) due to 
the deviations in the applied magnetic fields (SB = B — B , SB S = B s — B so ) 
can be expressed as: 

G(m t ,m s ,T,B,B s )=G°(m to ,m So ,T,B ,B so ) + AG. (7) 

Where G(m t , m s , T, B) is the free energy in the neighborhood of equilibrium, 
and G° is the equilibrium Gibbs free energy and AG is the production of the 
Gibbs energy due to the variance of the external field, and is given following 
form: 

AG = a(B - B ) + b(B s - B So ) + \c(m t - m to ) 2 + \d(m s - m SQ ) 2 
+\e(B - Bq) 2 + \f(B s - B SQ f + g(m t - m to )(m s - m so ) 

(8) 

+h(m t - m to )(B - B ) + i(m t - m to )(B s - B So ) 

+j(m s - m SQ )(B - B ) + k(m s - m SQ )(B s - B SQ ) 

The coefficients a — k are given in Appendix-B. 

In the sense of Onsager's theory of irreversible thermodynamics, the time 
derivatives of the antiferromagnetic and ferromagnetic order parameters are 
treated as generalized fluxes conjugate to their appropriate generalized forces. 
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One obtains the generalized forces (X mt ,X ms ) conjugate to the currents (rh t , rh s ) 
respectively, by differentiating AG with respect to (m t — m to ), (m s — m sQ ): 



dAG 

d(m t - m to ] 



X mt = — r = c(m t - m to ) + g(m s - m So ) + h{B - B ) + i(B s - B So ) , 



dAG 



d(m s — m 



so. 



Xm s = ~ — \ = d ( m * - m s ) + 9( m t - m to ) +j{H- H ) + k(B s - B, 

(9) 

The linear relations between the currents and forces may be written in terms 
of a matrix of phenomenological rate coefficients and since both m t and m s 
are odd variables under time inversion this matrix should be symmetric [35J: 



(10) 



Consequently, this matrix equation can be written in component form using 
equations ([9]), namely a set of two coupled, linear inhomogenous first order 
rate equations, 



m t 




Jmt 7 




X mt 






7 7m s 




Y 



m t = {^ mt c + jg)(m t -m to ) + {^ mt g + jd) (m s - m so ) 
Hlm t h + jj){B - B ) + ( 7mt i + ~fk)(B a - B S0 ), 
m s = (7c + "f ma g){m t -m t0 ) + (75 + "f ms d){m a - m so ) 



(11) 



+(7^ + lmJ)(B - B ) + (ji + i ms k){B s - B S0 ). 

In order to find the relaxation times, one considers the corresponding homo- 
geneous equations resulting when there is no external stimulation, namely 
B = B , B s = B So . Eqs. (11 ) then become 



'so- 



lm t c + ig im t 9 + ld 
ic + ~f ms g ^g + imj 



m s - m So 



(12) 



In iron group dihalides and all metamagnetic spin systems, there exists 
two order parameters, total and staggered magnetization (m t and m s ) which 
characterize the magnetic behavior of the system. Further, one can see from 
Eq.Q that, antiferromagnetic and ferromagnetic order are coupled to each 
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other. Due to this fact, there is an inteference between the ferromagnetic and 
antiferromagnetic relaxation processes. In the Onsager's theory of irreversible 
thermodynamics, the effect of cross effects between these two relaxation pro- 
cesses in the iron group dihalides is embedded in the kinetic rate coefficient 
7. Further, as it is discussed extensively by Barry and Harrington in Ref . [34] . 
other operator quantities should be added to the Ising Hamiltonian which do 
not commute with a z thereby introducing transitions within the spin system 
permitting, consequently, longitudinal relaxation. It is well known that, lon- 
gitudinal relaxation is closely related to spin-lattice relaxation time so that 
these added operator quantities should contain some kind of spin-lattice cou- 
pling, e.g. the frequently designated T\ appearing in the longitudinal Bloch 
equation representation of spin-lattice relaxation in solids |34j . 



The kinetic equations of motions given by Eqs.(ll) can be solved by as- 
suming a form for the solution ( m t — m to , m s — m So ~ ) by making use 
of the secular equation given below Thus, one obtains the following secular 
equation: 

\ + 7m« -76 1 C ~ Imb 

the resulting relaxation times can be found as: 



0. 



(13) 



1 1 , 1 l A l/2 

— = -19 ~ ^lm s d - ~7 mt c + -A 1 

T\ III 



1 1 , 1 I.,/, 

— = -19 ~ -7,lm s d - - 7mt c - -A 1 
T2 III 



(14) 



where A is 



A = (4 7ms d 7 # + 4 7 # 7mt c + ^d 2 - 2^ mt c^ ms d + ^c 2 + ^ mt g 2 lm s + ^l 2 dc 

(15) 

Temperature dependencies of the relaxation times for non zero external mag- 
netic field in the metamagnetic (antiferromagnetic) and paramagnetic phases 
for the cases {imlm s ~1 2 ~ 7m 7m J and (7 2 > 7 m 7mJ were discussed in Refs. 
12511261 . 



3.1 Derivation of Kinetic Equations Leading to AC Total Susceptibility 



If one stimulates the metamagnetic Ising system by magnetic field B os- 
cillating at an angular frequency u, all quantities will oscillate near the equi- 
librium state at this same angular frequency: 

m t -m t0 =m tl e I "\ m s - m SQ = m si e luit , B-B = B 1 e I " t . (16) 
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Substituting these equations into the kinetic equations given by Eqs.(ll) for 
B s = B SQ we find following form: 



'-Ioj + 7 mt c + jg) m tl + (j mt g + yd) m Sl = -y mt hB 1 , 



(7c + lm s g) m h + {-Iu + ig + im a d) m Sl = -yhB 1 



(17) 



After some algebra, we can find the matrix form of 



By 



-lm t h lm t g + ld 

m tl ~lh -Iui + yg + lm a d 

+ lm t C + 19 lm t 9 + id- 
lC + lm a g -/W + ig + ImJ 

The determinant in the denominator of Eq. (20) is the same as the secular 



determinant used to calculate the reciprocal relaxation times given by Eq.(13) 
except for the replacement of A by iu>. Hence Eq. (20) may be written as 



rati _ h (lmt{-Iu) + rf(7m t 7m s ~ T 2 )) 

Bi ~ 



-IUJ-T^ 1 ){-IUJ-T^) 



(19) 



Eq.(19) can be used in order to calculate the complex total magnetic suscep- 



tibility Xt(u)- This may be seen as follows: The induced magnetization (total 
induced magnetic moment per unit volume) of the spin-- Ising model with nn 
and nnn interactions is given by 



m t — m t 



Re(m tl e 



Iu>t\ 



(20) 



where mt x is the magnetization induced by a magnetic field oscillating at u. 
Also, by definition, the expression for Xt{<^) may be written 



m t - m to 



Re {xt{u)B x e 



Iujt 



(21) 



where 



Xt(u J ) = x' t (u J )-iXt(u J ) (22) 
here x m x'L ( u ) are the magnetic dispersion and absorption factor, re- 



spectively. Comparing Eqs. (21) and (22), one may write 

rrhi 



Xt{v) 



Si' 



(23) 
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and 



xiM = h- 



(tit 2 ) / 7 mt + Tid(7 



r 2 



7mt + T"2^(7 m s 7" 



J"2 - Ti , 



1 + W 2 Ti 2 



1 + w 2 r 9 2 



(24) 



(T1T2) / Ml (lm t + Tid(7 T 2 )) u;r 2 (7 mt + r 2 d(7 ma 7 mt - T 2 ))' 



(r 2 - 



1 + uHl 



1 + w 2 r 9 2 



(25) 



3.2 Derivation of Kinetic Equations Leading to Dynamical Staggered Mag- 
netic Susceptibility 



Now we assume that the system is stimulated by a staggered magnetic 
field B s oscillating at an angular frequency u, all quantities will oscillate near 
the equilibrium state at this angular frequency: 



m t - m to = m tl e Iu}t , m s - m so = m Sl e Iuit , 



B„ — B 



S ^SQ 



B Sl e 



Iu>t 



(26) 



embedding Eqs.(26) into the kinetic equations given by Eqs.(ll) for B = B 
one obtains 



(-Iu + 7 mt c + 7</) m tl + (ry mt g + jd) m si = -jkB Sl . 



(7c + lm s g) m tl + (-/w + 79 + 1mA) m si = -7m s kB Sl . 



solving the Eq(27 



yields the following matrix form of 



-Iu + 7 mt c + jg -jk 



7U 



B„ 



7C + lm a 9 



~lm s k 



m si 
B 

-iu + jmtC + jg im t g + i d 

lC + l ms g -IU + 70 + lm a d 

Following the same analysis given in Sec. 3.1 one obtains 

(lm s (-Iu) + C(7m t 7m s ~ T 2 )) 



m 
B, 



Sl 



-k- 



—Iuo — r x )(—Icu — t 2 



(27) 



(28) 



(29) 
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The induced staggered magnetization due to the effect of the oscillating stag- 
gered magnetic field is 



m q — m. 



Re(m Sl e 



(30) 



where m Soo is the staggered magnetization induced by a stagerred field oscil- 
lating at infinite frequency. Also, by definition the expression for the complex 
stagerred susceptibility reads 



m s — m Sc 



Re(xs(u)B sl e 



Iuit\ 



By making use of Eqs.(28) and (31) one obtains 



S 1 



B S1 ' 



(31) 



(32) 



finally, by using the relation Xs(w) = Xsi*- ) ~ %'( w ) we find 



(rir 2 ) / 7m s + ric(7 ms 7 mt - T 2 ) 7 ms + r 2 c(7 



r 2 



(r 2 - n) 



1 + u 2 rf 



1 + W 2 To 2 



(33) 

(nr 2 ) / wn(7 m + r!c(7 ma 7 mt ) - T 2 ) wr 2 (7 m + r 2 c(7 ms 7 mt - T 2 ) 



(r 2 - n) 



1 + w 2 ^ 2 



1 + W 2 T 9 2 



here x s ( u )y X s ( u ) are the magnetic dispersion and absorption factors respec- 
tively. 



4 Results 



Fig.2(a)-(b) shows the critical behavior of staggered magnetic dispersion, 
X s (w) and direct magnetic dispersion Xt( u ) with varying temperature in the 
low frequency region wr 2 <C 1 . One can observe from Fig. 2 (a) that real 
part of the staggered ac susceptibility (x s (u;)) rises rapidly with increasing 
temperature and tends to infinity near the second order phase transition. In 
accordance with the expectations the staggered dispersion factor converges to 
static staggered susceptibility for ur 2 <C 1 (compare Fig. 1(a) and Fig.2(a)). 
Fig2.(b) represents the temperature variation of the direct magnetic disper- 
sion factor in the neighborhood of T/v for B — 1.0 and rj = 2.0. The real part 
of direct complex susceptibility increases slowly with changing temperature 
and makes a cusp behavior near the second order phase transition. In addi- 
tion, Xt( u ) converges to static direct susceptibility in the low frequency region 
(compare Fig.2(b) and Fig. 1(b)). Moreover the behavior of Xs( u ) an d Xti u ) 
are independent of the frequency. This behavior is in accordance with results 
obtained for spin-1 Ising model [TT] . 
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Fig. 3(a)-(b) illustrates temperature variations of staggered magnetic dis- 
persion and direct magnetic dispersion for the high frequency region at several 
values of u in the neighborhood second phase transition point. has a 

local maximum before the second order phase transition point and shows a 
local minimum near the critical point in the high frequency region. In these 
figures the number accompanying each curve denotes the value of o and the 
vertical arrows refer to the critical temperature for B — 1.0 and rj = 2.0. 
The direct dispersion factor exhibits a maxima for T — > whose amplitude 
increases with rising frequency. This result is in parallel with findings for the 
total magnetic dispersion factor for the antiferromagnetic Ising model [3lj . 

In Fig.4(a)-(b) temperature variations of staggered magnetic absorption 
and direct magnetic absorption are given for the low frequency region at sev- 
eral values of u. x s exhibits a divergence, while Xt m &kes a local maximum 
near the critical point. Fig. 5 (a)- (b) shows temperature variations of staggered 
magnetic absorption and direct magnetic absorption for the high frequency re- 
gion at several values of u in the neighborhood second phase transition point. 
X s exhibits a local maximum at the phase transition point and amplitude of 
the maximum changes with frequency (See Fig. 5(a)). On the other hand, Xt 
makes local minimum at the phase transition point, but it exhibits a local 
maximum before the critical temperature (See Fig. 5(b)). The divergencies of 
X s (w) and Xsi 00 ) f° r WT 2 C 1 on both sides of the critical point are character- 
ized by the critical exponents. Thus, we may assume that, for temperatures 
smaller than the transition temperatures, x s ( u ) an d Xsi 00 ) follow laws of the 
form x s ( w ) ~ \T — T N \ lac and Xs( u ) ~ \T — T N \ lac , respectively. Here 7^, and 
7„ c are the dynamic critical exponents of the staggered magnetic dispersion 
and absorption factors. In order to calculate the mean field values of 7^, and 
7" c , we have sketched the the plots of x s ( u ) versus log(l — ?f-) and Xsi 10 ) 
versus log{l — ^) which are represented in Fig.6 (a). We have found only 

one linear part on these log — log plots for xl( w ) an d xl'( w )> i- e -, one values 
of 7^ c and 7^,, respectively. For each uj value the slope of the line is equal 
to —1 for j' ac ] and —2 for 7^,. These values are in accordance with the di- 
verging behavior of the magnetic dispersion and absorption factor near the 
second order phase transition temperature for B = Bq and B s = B So = 0. It 
is important to note that these findings are agreement with the critical ex- 
ponents found for the magnetic spin-1 Ising model [11]. Further these results 
are in well agreement with the MR studies of an AB type Ising model for for 
Ising antiferromagnetism [31] and Ising ferromagnet [SB]. In addition, Fig. 6(b) 
and (c) show logarithmic plot of staggered magnetic dispersion and staggered 
magnetic absorption versus reduced temperature (1 — T/Tc) for T < Tc at 
several values in the low frequency region and high frequency region. Finally, 
Table- 1 represents mean field critical exponents for staggered magnetic dis- 
persion, staggered magnetic absorbtion, direct magnetic dispersion and direct 
magnetic absorbtion in the low and high frequency regions. 
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5 Concluding Remarks 

Within the mean field approximation, we have analyzed steady state solu- 
tions of the spin- 1 metamagnetic Ising model under a time-dependent oscillat- 
ing external physical and staggered magnetic fields. In this paper, the formu- 
lation is based on a method which combines the equilibrium statistical theory 
of critical phenomena with the theory of irreversible thermodynamics. It is as- 
sumed that the amplitude of both the physical and staggered fields are so small 
that we made use of linear response theory in studying the magnetic relaxation 
processes in a metamagnetic system. In this system there exists two coupled 
relaxation processes which correspond to the relaxation of ferromagnetic and 
antiferromagnetic order parameters. We have shown that these processes are 
characterized by two distinct relaxation times (tj and T2 ). T2 is the dominant 
relaxation time which characterizes the critical slowing down of the staggered 
magnetization and therefore the temperature variation of the second relax- 
ation time determines the separation of the so-called low- and high-frequency 
regions. We should note that similar behavior has been observed in the in- 
vestigation of the antiferromagnetic Ising model with the same method [34]. 
Since — — > as T — > T N and keeping the frequency lo fixed, we observed 
the low-frequency behaviors followed by the high-frequency behaviors for the 
dynamic susceptibilities. One can see from Figs. 2(a) and 3(a) that the slope 
of the staggered magnetic dispersion curve chances in sign as T — > Tjy (posi- 
tive slope for the low-frequency region, negative slope for the high-frequency 
region). Similar behavior exists also in the high and low frequency regions 
for total magnetic dispersion factor. Finally, it should be emphasized that we 
have assumed in this study the rate coefficients have negligible temperature 
dependence. The validity of this assumption should be testified either by ex- 
periments or a more powerful theory such as path probability method. Kikuchi 
has represented an investigation of Order-disorder configurational relaxation 
on a bcc AB-type lattice and showed that The diagonal Onsager coefficients 
tend to finite values whereas the off-diagonal coefficient which characterizes 
the interference between coupled irreversible processes (37] tends to vanish as 
temperature approaches critical temperature [38] . 
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7 Appendix-A 

The list of coefficients in Eqs.(|6]): 

2cosh(2^)/ 2sinh(2^) 2 / 



an = 1 + 



k B T (cosh(2^) + cosh(2^)) kB T (cosh(2^) + cosh(2^ ^ 



2sinh(2^)sinh(2to^)t 



Qi2 — : > • 

A; B T (cosh(2 i y^) + cosh(2 i? g^) 

2s inh(2^±^)sinh(2^)/ 
a 2i — ; ~ ; \ 2 ' 



k B T (cosh(2^^) + cosh(2^ 



k B T 



(35) 

2cosh(2^±f^)i 2sinh(2^±f^) 2 t 



022=1- - -7 7-^77^ — ^—7 Y + 



fcflT ( cosh(2 ^/) + cosh (2^gH)) fc B T (cosh(2^) + cosh(2 



B s +m s f 
k B T ■ 



2 • 



2sinh(2^ii) s i n h(2^±^ 



k B T >_ ^""V k B T I 

c i — ; ~ : 1 ~r~2> 



k B T (cosh(2^^) + cosh(2 



B s +m a t 
k B T 



2cosh(2^gf^) 2sinh(2^±f^) 2 

C9 



k B T (cosh(2^) + cosh(2^)) kB T (cosh(2^) + cosh(2 



6 =1H 2cosh(2^g% 2sinh(2^)V 

fc B T(cosh(2^)+cosh(2tol)) (cosh(2^) + cosh(2^)) 2 A; B T' 



2sinh(2^^) sinh(2^g^)t 



cosh(2^f^) + cosh(2^±f^)) k B T 



>V2 



2sinh(2^)(sinh(2^^)/ 



(cosh(2^^) + cosh(2^±^)) fc B T 

(36) 



14 



2cosh(2^g^)t 



+ 



2s inh(2^g^) 2 t 



k B T (cosh(2^) + cosh(2^)) (cosh(2^) + cosh(2^)) 2 k B T 



2cosh(2^; 



2sinh(2 



B-m t f \2 
k B T > 



k B T (cosh(2^^) + cosh(2 



B s +m s t 
k B T 



k B T (cosh(2^^) + cosh(2 



B s +m s t 
k B T 



2 ■ 



2sinh(2^)sinh(2^) 
(cosh(2^) + cosH2^)f k B T 



(37) 
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8 Appendix-B 



The list of coefficients in Eq.(8): 

° = (f§) =-N»m 

\ / eq 



(d 2 G 

c = 



2 ' 



dm 



= X ( / - I - '/./-V 



4 V 1 + m s + m t 1 — m s — m t 1 + m t — m s 1 — m t + m. 



, /<9 2 G 
a = 



dm 2 . I 

s / eq 



-X{J + J' + \k B T 



4 V 1 + m s + m t 1 — m s — m t 1 + m t — m s 1 — m t + m s 



(d 2 G\ 

\ /eg 



; Ws) eq ' 



5 



\dm t dm s J 



-h B TN< 



4 V 1 + m s + m t 1 — m s — m t 1 + m t — m s 1 — m t + m s J ' 
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Fig. 1. Constant field (B) cross-sections of the staggered magnetic susceptibility (a) 
and direct susceptibility (b) in the neighborhood of critical point while B = 1.0 and 

'/ — z'J' ■ 




Fig. 2. Temperature variations of staggered magnetic dispersion (a) and direct mag- 
netic dispersion (b) for the low frequency region at several values of u in the neigh- 
borhood second phase transition point. 
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Fig. 3. Temperature variations of staggered magnetic dispersion (a) and direct mag- 
netic dispersion (b) for the high frequency region at several values of u in the 
neighborhood second phase transition point. The number accompanying each curve 
denotes the frequency. 




Fig. 4. Temperature variations of staggered magnetic absorption (a) and direct 
magnetic absorption (b) for the low frequency region at several values of uj in the 
neighborhood second phase transition point. 
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Fig. 5. Temperature variations of staggered magnetic absorption (a) and direct 
magnetic absorption (b) for the high frequency region at several values of uj in the 
neighborhood second phase transition point. The number accompanying each curve 
denotes the frequency. 
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1- 



-1- 

f-2- 

i-3 
-4- 
-5 



B=1.0 
11=2.0 



-8.5 -8.0 -7.5 -7.0 -6.5 -6.0 -5.5 
log(1-T/T) 



Fig. 6. Logarithmic plot of staggered magnetic dispersion and staggered magnetic 
absorption vs. reduced temperature (1 — T/Tq) for T < Tq at several values in the 
low frequency region (a) and high frequency region (b). Dotted lines are for x s an d 
solid lines for x s ■ 



23 



Table 1 

Mean field critical exponents for staggered magnetic dispersion, staggered magnetic 
absorbtion, direct magnetic dispersion and direct magnetic absorbtion in the low 
and high frequency regions. 



Quantity 


Frequency Range 


Exponent Value 


Critical Singularity 




2.1(T 7 - 1.10- 6 (low freq.) 


-1 


divergence 




2.10- 7 - 1.10~ 6 (low freq.) 


-2 


divergence 


xlM 


0.2 - 1.0 (high freq.) 





cusp 




0.2 - 1.0 (high freq.) 


+1 


cusp 




0.2 (high freq.) 


0.02 


cusp 


XtH 


0.4 - 1.0 (high freq.) 





cusp 




0.2 - 1.0 (high freq.) 





cusp 
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